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Abstract 


Cumulus clouds can become tilted or elongated in the presence of wind shear. 
Nevertheless, most studies of the interaction of cumulus clouds and radiation have as- 
sumed these clouds to be isotropic. This paper describes an investigation of the effect 
of fair-weather cumulus cloud field anisotropy on domain-averaged solar fluxes and 
atmospheric heating rate profiles. A stochastic field generation algorithm was used to 
produce twenty three-dimensional liquid water content fields based on the statistical 
properties of cloud scenes from a large eddy simulation. Progressively greater degrees 
of x-z plane tilting and horizontal stretching were imposed on each of these scenes, 
so that an ensemble of scenes was produced for each level of distortion. The result- 
ing scenes were used as input to a three-dimensional Monte Carlo radiative transfer 
model. Domain-average transmission, reflection, and absorption of broadband solar 
radiation were computed for each scene along with the average heating rate profile. 
Both tilt and horizontal stretching were found to significantly affect calculated fluxes, 
with the amount and sign of flux differences depending strongly on sun position rela- 
tive to cloud distortion geometry. The mechanisms by which anisotropy interacts with 
solar fluxes were investigated by comparisons to independent pixel approximation and 
tilted independent pixel approximation computations for the same scenes. Cumulus 
anisotropy was found to most strongly impact solar radiative transfer by changing the 
effective cloud fraction, i.e., the cloud fraction when the field is projected on a surface 
perpendicular to the direction of the incident solar beam. 
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1. Introduction 


The effect of three-dimensional (3D) cloud variability on the transport of solar flux in the 
atmosphere has been an active area of research for some three decades ( e.g ., McKee and Cox 1974; 
Aida 1977; Davies 1978; Harshvardhan and Thomas 1984; Welch and Wielicki 1984; Barker and 
Davies 1992; Cahalan et al. 1994b; Marshak et al. 1995b; Barker et al. 1998; O’Hirok and Gautier 
1998; Zuidema and Evans 1998; Benner and Evans 2001; Scheirer and Macke 2003; Di Giuseppe 
and Tompkins 2005). The 3D radiative transfer effect on domain- averaged solar fluxes has been 
divided into two physical processes, as summarized by Varnai and Davies (1999). The first, which 
is termed the “one-dimensional (1-D) heterogeneity effect,” arises from the nonlinear relationship 
between cloud optical depth and albedo. The mean transmission of a cloud with horizontally 
varying optical depth is more than the transmission of a uniform cloud with the mean optical 
depth. As a result, computed fluxes will be in error if a region with horizontally varying cloud 
properties is treated as homogeneous. The second process, which may be termed the “horizontal 
transport effect,” involves the horizontal propagation of radiation. Photons striking the top of a 
cloud can leak out the sides or diffuse out the top. For low solar elevations, photons can illuminate 
the side of a cloud, shading other parts of the field. Conversely, photons exiting cloud sides can 
enter the sides of other clouds and be scattered out their tops. All of these processes can affect a 
cloud field’s domain- averaged radiative properties, such as reflectance and transmittance. 

The independent pixel (or column) approximation (IPA) (Ronnholm et al. 1980; Cahalan 1989) 
was devised to treat cloud horizontal heterogeneity by integrating ID (plane-parallel) radiative 
transfer over the distribution of optical depths in a scene. The IPA has been found to be an accu- 
rate approximation for the computation of solar fluxes in overcast stratiform clouds (Cahalan et al. 
1994b; Chambers et al. 1997; Zuidema and Evans 1998). However, even if ID heterogeneity is 
accounted for, the horizontal transport effect causes potentially large solar flux errors in cumulus 
cloud fields (Barker 1996; Chambers et al. 1997; Benner and Evans 2001). The horizontal trans- 
port effect can be divided into two components. One component results from the way the apparent 
cloud geometry changes as a function of solar zenith angle. For example, the effective cloud frac- 
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tion (Harshvardhan and Thomas 1984) of cumulus clouds tends to increase with solar zenith angle 
because the sun illuminates the sides and tops of the clouds. The tilted independent pixel approxi- 
mation (TIPA) (Varnai and Davies 1999) attempts to account for this effect by performing radiative 
transfer on independent columns oriented along the slant path of incoming solar photons. This is 
equivalent to applying the IPA to a distribution of cloud optical depths adjusted for the angle of 
solar incidence. The second component of the horizontal transport effect is photon transport per- 
pendicular to the incident solar direction. This component includes radiative smoothing (Marshak 
et al. 1995a) and cloud side leakage (Varnai and Davies 1999). Perpendicular photon transport can 
only be correctly computed using full 3D radiative transfer techniques and thus requires knowledge 
of the complete structure of cloud fields. 

An ongoing problem in the field of 3D radiative transfer in the atmosphere has been the lack of 
three-dimensional cloud information. Most measurements are limited to the horizontal plane (e.g., 
satellite) or vertical cross sections (e.g., radar). The lack of 3D cloud structure data has limited 
our ability to accurately assess the real world importance of 3D radiative transport. It has also 
limited the types of cloud structures that have been investigated. One aspect of cloud geometry 
that has largely been ignored is anisotropy, i.e., structural variability that depends on the sampling 
direction. This oversight requires attention because anisotropy is common in broken cloud fields. 
Individual clouds may be elongated in a horizontal direction or tilted out of the vertical. On the 
larger scale, cloud fields frequently exhibit preferred directions, the most extreme example being 
cloud streets. 

Barker and Davies (1992) briefly addressed the importance of horizontal anisotropy in broken 
cloud fields to radiative transfer in climate models. Using a test cloud field based on a strongly 
anisotropic AVHRR scene, they found that reflectance varied by as much as 5% as a function of 
solar azimuth angle. They also detected smaller differences between the azimuthally averaged 
reflectances computed for the anisotropic field and for statistically similar isotropic scenes. They 
concluded that isotropic clouds may be a sufficient approximation for anisotropic fields in climate 
models. This study had two weaknesses. First, only one anisotropic scene was examined. Second, 
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the cloud fields on which these tests were run were essentially two-dimensional (only a single pixel 
thick) and had a uniform optical depth, r. Since more recent studies (Cahalan et al. 1994a; Titov 
1998; O’Hirok and Gautier 1998; Zuidema and Evans 1998; Varnai and Davies 1999) have shown 
that both horizontal r variations and vertical cloud structure are important to radiative transport 
in the cloudy atmosphere, further investigation of the radiative impact of horizontal anisotropy is 
required. 

Barker (1994) also examined the effect of shear-induced tilt on the albedo of cumulus clouds. 
This study was based on regular arrays of cuboidal clouds having a uniform volume extinction 
coefficient. The clouds were tilted to selected angles by consecutively offsetting the vertical layers 
in the horizontal direction. The albedos of the resulting cloud fields differed from the original 
field albedos by as much as 25%. These differences depended strongly on solar azimuthal angle 
as well as cloud tilt and were attributed primarily to changes in the effective cloud fraction. Since 
integration over all azimuthal angles reduced the total albedo difference to 5-10%, it was speculated 
that sheared cloud structures may not be significant in a climatological sense. This study relied 
primarily on model clouds with idealized geometric shapes and did not explore the effects of tilt 
on atmospheric heating rate profiles, leaving many questions unanswered. 

The current study attempts to improve our understanding of the impact of cumulus cloud field 
anisotropy on solar radiative transfer. In addition, we would like to elucidate the mechanisms 
through with these effects occur. Our ultimate purpose is to determine when, or whether, cloud 
field anisotropy must be accounted for. The input to this study is a set of stochastically generated 
cloud fields upon which varying degrees of anisotropy have been imposed. A large number of 
fully three-dimensional fair-weather continental cumulus scenes are considered to give statistical 
weight to the results. 
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2. Method 


a. Background 

1) Anisotropy measures 

In this paper, we consider two orthogonal types of anisotropy. The first, which we term hori- 
zontal anisotropy, corresponds to the “out of roundness” or elongation of cloud field structures in 
a horizontal plane. The second is tilt, or the degree to which the cloud structure shifts horizontally 
as a function of height. Separate measures are used to quantify these types of distortion. Since 
the methods used to compute these measures were presented in an earlier paper (Hinkelman et al. 
2005), we limit our discussion to the interpretation of these two variables. 

The anisotropy parameter (AP) is used to describe the degree to which the structures in a two- 
dimensional field vary as a function of direction. An individual cloud element is horizontally 
isotropic if its cross section in the horizontal plane is circular. It is anisotropic if the cross section 
is elliptical or linear. Anisotropy can also occur due to the arrangement of multiple cloud elements. 
For example, the elements could be spread randomly (an isotropic arrangement), linearly (highly 
anisotropic), or more closely packed in one direction than the other (somewhat anisotropic). The 
anisotropy parameter is sensitive to both of these types of “orientedness.” 

The anisotropy parameter consists of an amplitude and an angle. It can be calculated as a 
function of spatial scale or include a range of scales. In this work, it is the average of the anisotropy 
values for scales from 140 m to 1130 m (i.e., the smallest resolvable scale in the LES model to 
roughly 20% of the domain size). The significance of an AP value depends on the number of 
elements in the 2D field and can be evaluated using a statistical significance test as described in 
Hinkelman et al. (2005). AP amplitudes range from 0 for an isotropic field to 1 for a field that is 
perfectly linear. The AP angle indicates the direction of the long axis of the field’s structure. We 
specify AP angles according to meteorological conventions, i.e., 0°or 360° denotes a north-south 
(vertical) axis of anisotropy while 90° or 270° refers to an east-west (horizontal) orientation. (Note 
that anisotropy angles are unique only over a range of 180° , since they decribe alignment parallel 
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to an axis rather than pointing in a single direction.) 

Vertical tilt is defined as the horizontal shift of cloud field structure per unit height. This value 
is computed for every pair of neighboring horizonal slices through the cloud field, then averaged 
over the entire depth of the field. Tilt is expressed as a unitless quantity (in m in ~ 1 ) or in degrees 
away from the vertical. 

2) Choice of scale 

In this study, we have chosen to focus on cloud characteristics on the scale of 10-100 m. 
This is the limit to which cloud structure can be routinely measured at this time. Radar and lidar 
instruments, which provide some of the most detailed continuous measurements of cloud structure, 
are currently limited to vertical resolutions on the order of 10 m. Solar radiation also interacts with 
clouds on scales of tens of meters. For these reasons we employed cloud fields with a domain size 
of 6.4 km and 66.67 m sampling in the radiative transfer computations presented here. 

b. Stochastic cloud fields 

In order to isolate the effects of cumulus anisotropy on solar radiative transfer, we require a 
set of cloud fields for which the horizontal and vertical anisotropy vary over a desired range while 
the one- and two-point statistics of liquid water content remain fixed. A number of liquid water 
content fields with varying degrees of anisotropy and tilt were available from previous large-eddy 
simulations (LES) of radiatively-forced fair-weather cumulus over Oklahoma (Hinkelman et al. 
2005). However, the horizontal anisotropy and tilt varied independently of other cloud properties 
in these fields. We thus elected to study cloud scenes with prescribed anisotropy values produced 
by a modified version of the stochastic field generation algorithm of Evans and Wiscombe (2004). 
Multiple output cloud fields with matching liquid water content distributions and correlation struc- 
ture were produced for each set of anisotropy values. Averaging the results of the radiative transfer 
computations over these ensembles of similar fields allowed us to minimize the impact of variations 
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in cloud structure in individual scenes. 


The cloud field generation technique of Evans and Wiscombe (2004) is an extension of earlier 
methods of producing random two-dimensional fields with a known power spectrum under the 
assumptions of translational invariance and isotropy (e.g., Barker and Davies 1992; Evans 1993). 
The basic technique consists of filtering, or multiplying point-by-point in Fourier space, Gaussian 
random deviates by the desired amplitude spectrum and then Fourier transforming the output back 
to real space. The one-point statistics of the resulting fields remain Gaussian. If desired, the output 
array can be exponentiated or thresholded to produce a broken cloud field, but this changes the 
overall power spectrum. 

Extending this technique to include the vertical dimension is nontrivial because the assump- 
tions of isotropy and translational invariance do not hold in this direction, i.e., cloud structure does 
depend on height. Producing the desired vertical structure requires accounting for the correlation 
of liquid water content values at adjacent levels. Evans and Wiscombe (2004) generalized the 
Fourier filtering technique to the third dimension by using empirical orthogonal functions (EOFs) 
to decompose the vertical structure into independent (orthogonal) basis functions. (Note: Fourier 
series are the EOFs for translationally invariant fields.) 

The Evans and Wiscombe (2004) technique can be summarized as follows: First, the average 
cross correlation function is computed for all possible level pairs from the binary cloud masks of 
all of the input scenes. (Cloud masks are chosen for the spatial structure statistics because the 3D 
radiative effects in cumulus clouds are believed to be mostly a function of cloud morphology.) The 
resulting ’’binary correlation matrix” as a function of horizontal distance is then converted to an 
equivalent ’’Gaussian correlation matrix.” The frequency domain representation of the ’’Gaussian” 
correlation matrix is obtained via Fourier transform, and then an EOF decomposition is performed 
on each Fourier term (i.e., spatial frequency). This orthogonal function representation of the origi- 
nal correlation matrix is used as the basis of the stochastic cloud scenes. 

To generate a new output field from this information, a field of Gaussian random variables with 
the same dimensions as the input correlation matrix is first generated. The random field is then 
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multiplied point-by-point by the filtering amplitudes, which are the square roots of the eigenvalues 
from the EOF decomposition. The field is then transformed back into real space via consecutive 
EOF and Fourier transforms. Finally, the individual Gaussian values are mapped into the original 
liquid water content probability density function of the input scenes to obtain a broken cloud field. 
This process can be repeated using additional sets of Gaussian random deviates to produce any 
number of cloud scenes that are statistically similar to the population of cloud fields from which 
the horizontal power spectra and their vertical correlation properties were derived. 

The reason that an equivalent ’’Gaussian correlation matrix” is used in place of the binary cloud 
mask correlation matrix is that the final mapping of stochastic Gaussian values to the original liquid 
water content distribution is nonlinear, and therefore does not preserve the power spectra of the 
input. The equivalent ’’Gaussian correlation matrix” is selected so that the output fields created 
using this matrix will have the approximate correlation function of the binary mask fields after 
they are transformed to liquid water content values. 

Several modifications were made to the Evans/Wiscombe algorithm for this study. First, the 
use of input statistics from 3D FES liquid water content fields instead of 2D slices from upward- 
looking radar measurements was accomodated. This was accomplished by computing the mean 
3D binary correlation matrix for the FES fields and then reducing this array to two dimensions by 
integrating the horizontal dimension azimuthally. The correlation matrix created in this way has 
the same form as correlation matrices from 2D radar slices. However, since the original correlation 
matrix is computed in three dimensions, it is no longer necessary to convert the power spectra from 
2D fields to appropriate spectra for 3D fields as described in Appendix 1 of Evans and Wiscombe 
(2004). In addition, any anisotropy that may have existed in the input fields is eliminated when 
the third dimension is lost. We are then able to incorporate horizontal anisotropy and tilt into the 
output fields in a controlled manner. 

Horizontal anisotropy was added to stochastic cloud fields by elongating the liquid water con- 
tent power spectrum for each vertical level in the North-South direction and compressing it in the 
East- West direction. This was done at the beginning of the cloud generation process, when the 
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EOF vectors were multiplied by the Gaussian random deviates in the two-dimensional array. At 
every point in the array, x, y, instead of multiplying by the EOF vector corresponding to radial 
spatial frequency bin 

k r = ^kl + k%, (1) 

where k x and k y are the bin indices corresponding to positions x and y, the EOF vector for spatial 
frequency bin k' r was used. Here 


k' = 


\ 


(. k x stretch ) 2 + 


( k y 

\ stretch ) 


( 2 ) 


and stretch is the parameter used to control the horizontal anisotropy. This produces a cloud field 
elongated in the East-West direction after the matrix is inverse Fourier transformed back to the 
spatial domain. 

Controlling the degree of horizontal anisotropy produced in the output fields is complicated 
by the fact that conversion from the Gaussian distribution back to the original LWC values is a 
nonlinear operation. As a result, the final output clouds will not have the same degree of horizontal 
anisotropy as was set by the input parameters. Repeated trials are thus required to determine the 
parameter values required to obtain the desired levels of anisotropy. 

Arbitrary vertical tilt from West to East was added by multiplying the Fourier amplitudes by a 
linear phase shift given by 

de = e^ /Cxdx)ni , (3) 


where c x is one half the total number of FFT frequencies, k x is the horizontal (W-E) spatial fre- 
quency index, and dx = iz tilt where iz is the index of the vertical cloud level and tilt is a 
parameter specified in horizontal grid spacings per vertical grid cell. 

Two three-dimensional liquid water content fields from each of six runs of the LES base case 
in Hinkelman et al. (2005) were used as input to the stochastic cloud field generation model. The 
six runs differed only in the random potential temperature fields used to initiate the convection. 
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These well developed fair-weather cumulus fields were all taken from simulation times between 
8 and 9 hours to provide a statistically similar ensemble of scenes. Vertical levels that were free 
of liquid water in all of the scenes were not included in the input to the stochastic cloud field 
program. For this reason, the output fields consisted of 96x96x37 arrays at the LES grid spacing 
of 66.66 x 66.66 x40.0 m 3 . Examples of LES scenes used as input to the algorithm are shown in 
Figure 1. 

Twenty realizations of isotropic cumulus cloud fields were generated using the stochastic 
method described above. Ten variations of each of these basic scenes were also produced by 
incorporating tilt or anisotropy. Five of the variations included tilts of 0.33, 0.66, 1.00, 1.33, and 
1.66 m m' 1 (corresponding to 0° , 18° , 33° , 45° , 53° , and 59° from vertical), while the other five 
had anisotropy parameters of approximately 0. 15, 0.25, 0.35, 0.45, and 0.55 over the middle spatial 
frequencies. Note that these values were computed for individual scenes rather than as averages 
over 10 scenes as in Hinkelman et al. (2005). The anisotropy values used here were selected to 
cover a range from zero to ~ 1 60% of the maximum single scene value in the data from Hinkelman 
et al. (2005). Scenes with both tilt and horizontal anisotropy were not created so that it would be 
possible to isolate the effects of these two types of distortion. The degree to which we were able 
to control the anisotropy of the stochastic cloud fields is shown in the graphs in Figure 2. 

A sample sequence of stochastic cloud fields with increasing levels of imposed horizontal an- 
isotropy is shown in Figure 3. The corresponding anisotropy parameter is indicated above each 
field. The anisotropy parameter value of 0.08 indicates the ensemble of nominally isotropic fields. 
This value is nonzero because some anisotropy occurs by chance in natural-looking cloud fields. 
Prior LES studies indicate that anisotropy parameter values of cloud fields typically range up to 
about 0.4. Higher anisotropy parameters would occur only in extremely oriented fields, such as 
when cloud streets are present. The sequence of tilted clouds made using the same random com- 
ponents as the anisotropic fields of Figure 3 is shown in Figure 4. Note that imposing tilt in the 
West-East direction leaves the North-South direction unchanged. Mean (standard deviation) cloud 
properties for the isotropic scenes are cloud fraction of 0.24 (0.05) and cloud optical depth of 19.2 
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(5.0). All of the clouds are less than 1450 m thick. 

A comparison of the average cloud fraction as a function of height for the ensembles of stochas- 
tic fields with differing degrees of tilt and anisotropy is presented in the upper portion of Figure 5. 
The ensemble-averaged cloud fraction varies so little as anisotropy increases that all of the mean 
value curves are coincident. At the cloudiest heights, the standard deviations of cloud fraction are 
higher for the stochastic ensembles than for the original group of scenes. However, the fact that 
the standard deviation curves are identical for all of the sets of anisotropic scenes suggests that 
this increase is a result of the limited number of stochastic scenes employed. Similar behavior is 
seen in the statistics of average liquid water content per cloud layer shown in the lower portion of 
Figure 5. The average liquid water content at each height is also conserved as the tilt or horizontal 
anisotropy increases. The fact that the bulk cloud properties are so well conserved means that we 
will be able to make meaningful comparisons of the radiative fluxes and heating rates calculated 
for the fields as tilt and stretch are varied. 

c. Radiative transfer computations 

Monte Carlo radiative transfer calculations were performed for the twenty sequences of tilted 
or stretched cumulus scenes described above. The atmospheric property arrays required for the 
radiative transfer computations were created by inserting the LWC arrays from the stochastic cloud 
model into larger arrays defining the computational domain. The domain extended from the surface 
to 40.0 km with the cloudy layers between 1060 m and 2500 m. The cloud droplet effective radius 
in each cell was calculated from the local liquid water content value assuming a constant drop 
number density of 300 cm" 3 and a lognormal droplet distribution with a logarithmic width of 0.38 
(Miles et al. 2000). A single atmospheric profile of temperature and pressure was used throughout 
the entire domain. Average values from the base LES cases were applied up to about 4000 m. 
Above 4300 m, the McClatchey midlatitude summer atmospheric profile was used. In between, 
the values were interpolated to splice the two curves together smoothly. 

Solar radiative transfer through the cloud scenes was computed using a Monte Carlo model 
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(Benner and Evans 2001) based on the maximal cross section method of Marchuk (Marchuk et al. 
1980) and tested by comparison to other algorithms in ICRCCM-3 (Barker et al. 2003). The Monte 
Carlo algorithm included molecular Rayleigh scattering and incorporated molecular absorption 
via the shortwave rapid radiative transfer model (RRTM) correlated-k distribution (Mlawer et al. 
1997), comprised of 14 bands from 0.2 to 12.2 fim. Optical properties averaged over each wave- 
length band were calculated explicitly using Mie theory for the droplet distributions. Up to 2000 
Legendre coefficients were used to represent the phase functions for each of 40 effective radius 
values ranging from 0.5 to 20 fim in 0.5 //m intervals. The surface was chosen to be perfectly 
absorbing. 

Full three-dimensional Monte Carlo computations were made for each of the 220 scenes for 
three solar zenith angles (0° , 30° , and 60° ) and azimuthal angles of 0° and 90° . The azimuthal 
angles are defined in Figure 6. The 0° azimuth is defined as solar input from the east, which is 
parallel to the stretching and tilting axes with the tops of the tilted clouds pointing towards this 
direction. The 90° azimuth is broadside to the stretched and tilted clouds. A 270° solar azimuth 
angle was not employed, as this gives results equivalent to the 90° case for both the tilted and 
stretched fields. The 180° azimuth indicates solar input from the west. As for 0° , this orientation 
is parallel to the axis of stretching and tilting and yields similar radiative results for the stretched 
fields. However, the tops of the tilted cloud fields point away from this direction, so calculations 
with a 180° solar azimuth angle were performed for the vertically tilted cloud fields. Runs in 
which the photons were injected at random solar azimuth angles were also performed in order to 
investigate possible averaging effects due to variations in cloud-sun orientation over time or space. 

In addition to the full three-dimensional radiative transfer computations, independent pixel 
approximation (IPA) and tilted independent pixel approximation (TIPA) calculations were carried 
out for the same scenes and solar zenith and azimuth angle combinations. These calculations 
were performed using the Monte Carlo algorithm by neglecting photon translation perpendicular 
to the vertical (IPA) or sun-oriented (TIPA) atmospheric columns. The results of these approximate 
radiative transfer methods were compared to the output of the full 3D Monte Carlo to elucidate the 
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ways in which anisotropic fields affect solar radiative propagation. 

The variables provided by the Monte Carlo calculations were the domain-averaged downward 
direct and total surface fluxes, upward flux at the top of atmosphere, and total absorption. The 
diffuse ratio was computed as the difference between the total and direct surface fluxes normalized 
by the total downwelling surface flux. In addition, domain-averaged heating rate profiles were 
calculated. These heating rates were based on the average absorbed flux, temperature, and pressure 
at each model height. Six million photons were used for each Monte Carlo run. This number was 
selected based on accuracy and run time considerations. Tests run with a typical cloud field and a 
solar zenith angle of 30° yielded the error estimates listed in Table 1. 

3. Results 

a. Domain-average fluxes 

Domain- averaged transmission, reflection, absorption, and diffuse flux ratio as a function of 
cloud field horizontal anisotropy are presented in Figure 7. Each point on the curves represents an 
average over the appropriate ensemble of 20 fields, while the error bars indicate the standard error 
of the mean for each ensemble. The top left panel of Figure 7 shows that cloud field anisotropy 
does affect the transmission of solar radiation for nonzero solar zenith angles. This effect depends 
on the orientation of the incident beam relative to the axis of cloud field stretching. For solar 
input parallel to the long axis of the clouds the transmission increases, while it decreases for solar 
incidence broadside to the stretched fields. These effects are consistent with a change in effective 
cloud fraction due to the stretching for nonzero solar zenith angles. Stretching along the east-west 
axis increases the cloud cross section presented to solar fluxes traveling perpendicular to the x-axis 
and thus reduces the proportion of the solar flux reaching the surface. The opposite is true for 
fluxes with angles of incidence parallel to the stretching axis. Although it is not obvious from the 
graphs, the values in Table 2 demonstrate that the impact of the anisotropy is small but nonzero 
for random solar azimuth angles. This is because the increase in transmission for incident angles 
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parallel to the stretching axis is smaller than the decrease in transmission that occurs for broadside 
incidence. The fact that increasing horizontal anisotropy has relatively little effect for overhead 
sun confirms that cloud properties are well conserved as the anisotropy is increased. 

Results for reflection or albedo, shown in the upper right panel of Figure 7, are essentially 
the inverse of the transmission results. The albedo increases with increasing anisotropy when the 
solar flux is incident on the side of the clouds and decreases with increasing anisotropy when the 
incoming flux is parallel to the stretching axis. Changes in the total absorption are small (note 
the change in vertical scale and the values listed in Table 3), but vary in direction according to 
the effective cloud fraction, i.e., absorption increases/decreases with anisotropy for side/parallel 
solar incidence. The largest relative changes with anisotropy in Figure 7 occur in the diffuse ratio 
(bottom right panel). As we would expect, the diffuse ratio increases with increasing effective 
cloud fraction and vice versa. 

Domain- averaged transmission, reflection, absorption, and diffuse ratio as a function of cloud 
tilt are presented in Figure 8 and the transmitted flux values are listed in Table 4. In general, the 
effect of tilt is larger than that of stretching for most solar geometries. Even for a solar zenith 
angle of 0° , tilt significantly increases the albedo and decreases transmission. This is because 
shifting cloud layers horizontally causes an increase in cloud fraction for overhead sun, while cloud 
fraction is conserved in the stretching process. Apparent cloud fraction also increases with tilt for 
the 90° azimuth (broadside to the tilting axis). Thus the transmission decreases with increasing tilt 
from this direction as well. 

It is noteworthy that the transmission values calculated for solar incidence parallel to the 
stretching axis (azimuths of 0° and 180° ) are very different in the presence of tilted fields. For 
a 180° azimuthal angle, the clouds tilt away from the incoming radiation, thus presenting a greater 
cross section to the incident flux and decreasing solar transmission. In contrast, for the 0° azimuth, 
the clouds tilt towards the sun, which reduces their effective cross section. The maximum trans- 
mission for tilted cloud fields occurs at the 0° azimuth when the tilt of the clouds matches the solar 
zenith angle, so that the effective cloud fraction is minimized. For the 30° solar zenith angle, this 
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occurs at a tilt of about 0.66. The transmission peaks at this tilt increment, then decreases as the tilt 
continues to grow. No peak is observed in the 60° solar zenith angle transmission curve because 
the imposed tilts do not exceed 60° . 

In contrast to horizontal stretching, the impact of tilt on absolute fluxes generally diminishes for 
larger solar zenith angles. (See Table 4.) This is because the 60° solar zenith angle is nearly a side 
view; lateral displacements of cloud structures don’t change the total cross-sectional area much 
from this perspective. Only for the tilt-aligned 0° azimuth is the largest effect on the downward 
flux observed for the 60° solar zenith angle. 

As for the horizontally anisotropic fields, the changes in the albedo with tilt are nearly opposite 
the differences found for transmission. Again, the variations in total absorption, listed in Table 5, 
are small but nonzero. Except for the 0° azimuth, the diffuse ratio always increases with tilt, 
reaching values over 0.4 for the most tilted scenes at low sun angles. 

Again as with the stretched cases, use of random azimuthal angles does not “average out” the 
effects of tilt. The reason for this is evident from the differences in transmitted flux relative to the 
untilted scenes listed in Table 4. Increasing tilt decreases the transmission (increases the albedo) 
for all azimuthal angles except 0° and the effect at 0° is not enough to outweigh the contributions 
from the other solar azimuthal angles. 

From Tables 2 and 4 it is apparent that not taking horizontal anisotropy into account can cause 
errors in the downward flux at the surface close to 10 Wm~ 2 , while neglecting vertical tilt can lead 
to errors as high as 40 Wm~ 2 in extreme cases. Errors of this magnitude can be important to energy 
budget calculations if they persist for any length of time. The errors in the absorbed flux are much 
smaller in magnitude, but still significant as a fraction of the total fluxes. 

b. Atmospheric heating rate profiles 

Although the changes observed in domain-averaged absorption due to cumulus tilt or anisot- 
ropy were small, their spatial distribution could have important implications for the local energy 
budget or the dynamics of the cloud fields. For this reason, we analyze the heating rate changes 
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caused by cloud field anisotropy as a function of height in the domain. 

The change in heating rate profiles for cloud scenes as vertical tilt increases is shown in Fig- 
ure 9. A separate plot is shown for each solar zenith angle and azimuth angle combination. The 
heating rate profile for the ensemble of untilted case is illustrated by the black line. Differences 
between this profile and those for the progressively more tilted cases are shown in various shades 
of grey. The general shape of the base heating rate curve is determined by the distribution of liquid 
water with height. Since most of the absorption occurs within the first few optical depths of ab- 
sorbers encountered by the photons, the heating rate maximum occurs near cloud top. The heating 
rate is then fairly constant throughout the depth of the clouds before falling off rapidly below cloud 
base. 

From the perspective of atmospheric heating rates, the effect of tilt is to alter the apparent cloud 
top area. For example, for overhead sun, skewing enlarges the upper surface area of the cloud field. 
This causes heating rates to rise at all but the top levels because new upward surfaces are exposed at 
lower elevations as the sides of the clouds peak out from below. This same effect occurs for all solar 
azumithal and zenith angles except the azimuth of 0° . For this geometry, increasing tilt decreases 
the apparent cloud fraction, reducing the cloud top area and heating rates correspondingly. The 
occurrence of a position of minimum cloud fraction is apparent in the curves for a solar zenith 
angle 30° and azimuth of 0° . Starting from no tilt, the heating rates first decrease and then increase 
after the minimum in apparent cloud top surface area is reached at a tilt of 0.66. Although a similar 
decrease in heating rate with increased tilt is observed for the 60° solar zenith angle, a minimum 
is not reached due to insufficient tilt in the cloud fields. The maximum local heating rate change 
caused by tilt in this study is +0.65 Kday -1 , corresponding to an increase of 17%. This occurs at 
the top of the shoulder region of the profile (at 2.22 km) for 30° solar zenith angle and 180° solar 
azimuth angle. 

Whether the effects of tilt on atmospheric heating rate profiles can be cancelled out by varia- 
tions in azimuthal angle depends on the solar zenith angle. For solar zenith angles of 0° and 30° , 
heating rates increase with tilt as the apparent cloud fraction increases even for random azimuthal 
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angles. But at 60° solar zenith angle, there is a great deal of side incidence at all azimuthal an- 
gles, so the changes in heating rate with cloud tilt are of lesser magnitude. The strong decrease in 
heating rates with tilt for 0° azimuth more than offset the increase for 180° . With the addition of 
small heating rate increases for azimuths of perpendicular solar incidence, there is nearly no net 
change in the atmospheric heating rate profile at a 60° zenith angle when random azimuthal angles 
are employed. 

The effect of horizontal anisotropy on atmospheric heating rate profiles is illustrated in Fig- 
ure 10. Again, there is a separate graph for each zenith-azimuth angle combination and the shade 
of the lines indicates the level of anisotropy in the ensemble represented. Note that there are fewer 
graphs because the 180° azimuth is not considered. 

The behavior of the heating rate profiles with increasing horizontal anisotropy is as follows: 
For overhead sun, almost no change is observed because the cloud fraction and vertical optical 
depth distributions are well conserved as the clouds are stretched. For more oblique angles of 
incidence, heating rates decline for solar azimuthal angles parallel to the stretching axis (0° ). 
For these relative solar positions, the apparent optical depths rise but the cloud cross sections 
decline, giving less effective cloud top area. The opposite occurs for solar incidence perpendicular 
to the stretching axis where the apparent optical depths decrease but the cloud top surface area 
increases, yielding higher heating rates. For this reason, the changes in the heating rate profiles 
due to horizontal anisotropy tend to cancel for random azimuthal angles. Horizontal anisotropy 
has a significantly smaller impact on heating rate profiles than vertical tilt does for the scenes used 
in this study. The maximum heating rate profile change caused by horizontal anisotropy is only 
+0.25 Kday -1 . However, this occurs for a solar zenith angle of 60° and azimuth of 90° , so that the 
relative change is more than 10%. This maximum also occurs at 2.22 km. 
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4. Discussion 


a. Mechanisms by which anisotropy affects solar radiative transfer 

Comparing results for the independent pixel approximation (IPA), tilted independent pixel ap- 
proximation (TIPA), and full three-dimensional radiative transfer calculations illuminates the pro- 
cesses by which anisotropy affects solar radiative transport since each successive approximation 
further restricts the accuracy of photon tracing through the cloud scenes. Differences between the 
IPA and 3D results reflect both parts of the 3D effect - constraining the photons to stay within a 
given vertical column eliminates scattering into a neighboring column and neglects the changes 
in apparent cloud geometry for nonvertical solar illumination. TIPA, on the other hand, accounts 
for the apparent changes in cloud geometry as a function of solar incidence. Differences between 
TIPA and 3D results thus stem entirely from the neglect of scattering perpendicular to these paths. 

Differences between the transmitted fluxes for IPA and 3D computations for the stretched and 
tilted scenes are shown in the upper panels of Figure 11. Neglect of both types of “3D effects” 
leads to errors as large as 23 Wm -2 for the stretched cases and 29 Win 2 for the tilted cases. Since 
the 3D results are highly dependent on solar zenith and azimuth angles while the IPA is nearly 
insensitive to solar position, the IPA errors depend strongly on the relative sun angles. For the 
stretched scenes, the IPA errors for 0° and 90° tend to be offsetting, but even for random azimuthal 
angles, the IPA results still differ from the full 3D transmitted fluxes. The same is true for the 
0° and 180° azimuth results for the tilted fluxes. For both the tilted and stretched scenes, only the 
fluxes for 30° solar zenith angle and random (also 90° for the tilted scenes) azimuth angles stay 
close to 0 Wm~ 2 as the degree of anisotropy increases. 

The differences between the TIPA and full 3D transmitted fluxes, illustrated in the lower panels 
of Figure 11, indicate that changes in cloud geometry are the chief driver of the variations in 
domain- averaged solar radiative quantities that we have observed with increasing anisotropy. Since 
both the TIPA and 3D computations use the proper cloud geometry for each direction of solar 
incidence, the differences between their results are due primarily to omission of scattering out of 
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this direction. Figure 1 1 shows that the magnitude of the perpendicular transport is nearly constant 
regardless of tilt, horizontal anisotropy, and even relative solar position, having values between 9 
and 12 Wm~ 2 in all cases. This suggests that the perpendicular transport is largely a function of the 
basic properties (such as LWC, cloud amount, and size) of clouds in a scene and only minimally 
impacted by minor changes in cloud geometry. 

Changes in apparent cloud geometry due to relative solar position can be further divided into 
changes in effective cloud fraction and mean cloud optical depth. Cloud fraction is a measure of 
the extent of a cloud field while optical depth is related to the intensity or duration of the interaction 
between the solar photons and cloud droplets. Effective cloud fraction and effective optical depth 
were computed for each stochastic cloud field by Monte Carlo sampling of photons injected at 
the desired solar zenith and azimuth angles. For the correct (3D) relative position of the sun, the 
effective cloud fraction curves (not shown) precisely mirror the albedo curves of Figures 7 and 8 
while the mean cloud optical depths change in the opposite direction as anisotropy increases. Note 
that the observed inverse relationship between apparent cloud fraction and optical depth is a special 
characteristic of the fields used in this study that arises because the total liquid water content of the 
cloud fields is held fixed as anisotropy increases. 

Earlier, we argued that the changes in solar flux terms computed for cloud scenes with increas- 
ing levels of anisotropy are consistent with changes in effective cloud fraction. This hypothesis is 
supported by the plots of transmission and diffuse flux ratio as a function of effective cloud fraction 
and optical depth shown in Figure 12. Other things being equal, we would expect transmission to 
decrease as cloud fraction or optical depth increase. While this is the case for the effective cloud 
fraction, transmission is observed to increase as optical depth increases. This implies that, at least 
for the cloud scenes studied here, effective cloud fraction has a greater impact on transmitted fluxes 
than does cloud optical depth. For the diffuse flux ratio, the expected trends are opposite those for 
transmission, that is, diffuse flux ratio should increase as either effective cloud fraction or optical 
depth increase. Again, we find that the expectation holds true for cloud fraction but not for optical 
depth. By physical reasoning, then, the changes in mean cloud optical depth caused by increasing 
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cloud anisotropy cannot cause the observed changes in the solar fluxes. [Note that only results for 
the tilted scenes are shown in this figure. The trends for the stretched scenes are the same, but the 
cloud fractions and optical depths fall within a more limited range.] 

Although the effective optical depth seems to play little part in the changes in the domain- 
averaged radiative transfer associated with increasing anisotropy, it is relevant to the absolute 
magnitude of the domain-averaged fluxes as a function of solar zenith angle. The large separa- 
tion between the transmission curves as a function of effective cloud fraction in Figure 12 are 
explained by the overall increase in cloud optical depth as solar zenith angle increases. In contrast, 
the computed diffuse flux ratios have nearly a one-to-one correspondence with the corresponding 
effective cloud fraction values, suggesting that the diffuse flux ratios for the scenes considered here 
are almost entirely a function of effective cloud fraction alone. 

We conclude that tilt and horizontal anisotropy affect solar radiative transport primarily by 
changing the effective cloud fraction of a scene. Changes in the mean optical path are of lesser 
importance. Transport perpendicular to the direction of incidence is of a significant magnitude for 
the cases examined here, but is largely unchanged by simple rearrangements of the liquid water 
content distributions that leave the basic cloud properties unchanged. 

b. Potential approximate methods of radiative transfer for anisotropic cumulus clouds 

Given the results of the various radiative transfer computations described above, we can con- 
sider the effectiveness of various treatments of cumulus anisotropy in modeling and measurement 
situations. There are two factors that control the accuracy of such treatments: the amount of cloud 
structure information that is available and the degree of approximation in the radiative tranfer al- 
gorithm employed. If the complete 3D structure of the cloud field is known, IPA, TIPA, and full 
3D radiative transfer methods are all applicable. Of course the 3D computations will give accu- 
rate flux and heating rate values. However, if a faster method is desired, applying TIPA will be 
reasonably accurate. As discussed previously, the the error in transmitted fluxes computed using 
TIPA is always negative and on the order of 10 Wm -2 for the type of scenes used here. For other 


20 



cloud fields, the size of the bias could be estimated from known properties of the cloud fields. This 
would require previous computations on sample scenes with various values of mean cloud size 
and optical depth, cloud thickness, and average spacing between clouds. There is no advantage to 
using IPA because it gives widely varying and frequently large errors (-30 to 25 W m 2 for the 
scenes used here) while being no more computationally efficient than TIPA. 

In most modeling and even measurement situations, the full 3D cloud field structure is not 
known at high resolution. The basic information needed as input to TIPA for anisotropic cumu- 
lus fields is the optical depth distribution as a function of solar zenith and azimuth angles, i.e., 
P(t\ 9 0 , (f> 0 ). To be useful for climate modeling, how this optical depth distribution depends on 
large-scale thermodynamic and dynamic variables (e.g., water vapor concentration, convective sta- 
bility, and wind shear) would have to be known. Information about the P(r\9 0 , 0 O ) distribution 
could potentially be measured with instruments such as scanning microwave radiometers. A more 
likely approach for climate modeling would be to 1) express P(r\9 0 ,(j ) 0 ) in terms of relatively 
few parameters (e.g., by using gamma distributions in r), 2) run cloud dynamical models for a 
large number of forcing conditions, and 3) use the resulting model cloud fields to determine the 
parameters of the anisotropic optical depth distribution as a function of large-scale variables. (For 
example, see Hinkelman et al. (2005) for sample relations between wind speed, shear, and anisot- 
ropy.) 

c. Importance of anisotropy effects to various applications 

The results of this study indicate that not accounting for horizontal anisotropy or vertical tilt in 
cumulus cloud fields can cause errors of up to 40 Wm~ 2 in instantaneous downwelling solar surface 
fluxes at spatial scales on the order of 10 km. An important question is then whether anisotropy 
can significantly affect solar radiative parameters on longer temporal or larger spatial scales. In 
order for this to be the case, the following conditions would have to hold. First, the cloud fraction, 
optical depth, and aspect ratio of the cumulus clouds would need to be large enough for the field 
to produce significant radiative effects regardless of the anisotropy. In addition, the anisotropic 
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cloud field would have to exist over most of the period and region, the direction of the anisotropy 
would have to be uniform over the same intervals, and the relative orientation of the sun to the 
cloud features over the period of interest would need to not average out the relevant anisotropy 
effects. Based on the lifetime and extent of cumulus fields and midlatitude dynamic features, 
we might expect anisotropic fields to exist for at most 3-4 hours over a region on the order of 
100 km. It is also conceivable that the type and direction of anisotropy in such fields would be 
reasonably uniform. Given this, the important question becomes how much any temporal and 
spatial averaging would reduce the anisotropy’s radiative impact. Based on our results for random 
azimuthal angles, averaging would tend to minimize the effects of stretched clouds. However for 
tilted clouds, substantial radiative effects could remain. This is because the effective cloud fraction 
increases with tilt for all solar azimuthal angles except 0° . Thus cumulus tilt could noticeably 
change solar transmission and albedo along with the diffuse ratio for scales as large as 100 km 
and a few hours. In cases with more extensive averaging in time or space, anisotropy effects are 
probably insignificant. 

We can estimate the magnitude of the effect of anisotropy on transmitted solar flux for realistic 
cumulus fields over a longer period from the anisotropy values computed for LES cumulus scenes 
in Hinkelman et al. (2005). Table 5 in that paper shows anisotropy values for LES runs using a 
variety of typical environmental conditions vector averaged over 2.5 hr. The maximum anisotropy 
parameter value listed is 0.229, which occurred for a vertical wind shear of 4 m s -1 km -1 . Table 2 
indicates a maximum change in the downward solar flux of about -5 Wm~ 2 for a solar zenith 
angle of 60° at this level of stretching. The maximum tilt listed is 0.731, for the case with a 
uniform wind speed of 20 m s -1 . This tilt could decrease the downward solar flux by as much 
as 13.4 Wm~ 2 , again for a 60° solar zenith angle. Thus anisotropy could produce 5-12 Wm~ 2 
variations in the downwelling surface solar flux. The importance of such variations would depend 
on the application. 

Changes in atmospheric heating rate profiles due to tilt or horizontal anisotropy could alter the 
structure or microphysics of the clouds themselves. However, the impact of anisotropy on the de- 
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velopment of cumulus clouds is difficult to estimate from the analysis presented here. The heating 
rates computed here, which vary less than 1 Kday -1 as cloud structure changes, are averages over 
the entire domain, including the clear areas. A three-dimensional analysis of heating rate changes, 
which would identify “hot spots” (Vogelmann et al. 2001) within the clouds, would provide more 
insight into this matter. Large-eddy simulations with interactive 3D radiative transfer could also 
shed light on this subject. Given that the wet adiabatic cooling rate is about 50 K h '(Cotton and 
Anthes 1989), anisotropy would have to cause pockets of very strong heating to have any impact 
on cloud development. It is more likely that anisotropy-induced solar heating would be effective 
after cloud development, tending to shorten cloud lifetimes. Further study is required to evaluate 
these effects. 

d. Caveats 

It should be noted that the results described in this paper were obtained for a group of continen- 
tal fair-weather cumulus-type fields having similar properties and created by a stochastic model. 
All of the fields had similar cloud fractions and aspect ratios, and cloud depths did not exceed 2 km. 
However, an ensemble of twenty fields was used for each tilt or anisotropy state investigated. Our 
confidence in the reliability of the computed flux and heating rate trends for the type of cloud fields 
employed is therefore high. Nevertheless, the radiative effects of anisotropy and tilt in other types 
of cloud fields may differ. For example, we expect that cloud field distortions in cumulus with 
greater vertical development will produce larger changes in surface radiative fluxes while anisot- 
ropy effects for more overcast scenes are anticipated to be less significant. As in most studies of 
radiative transfer in the cloudy atmosphere, assumptions were made regarding cloud microphysical 
properties, specifically the number concentration of the droplets. Despite any infidelity due to the 
artificial nature of the test scenes, we expect the overall trends noted in this work to hold because 
the greatest effect of anisotropy on radiative processes was due to gross changes in the apparent 
shape of cloud boundaries. 
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5. Conclusions 


The radiative impact of horizontal elongation and vertical tilt of fair-weather cumulus clouds 
has been assessed. Statistics describing the structure of cumulus scenes produced by a large eddy 
simulation model were used as input to a stochastic field generation algorithm, which created a 
series of 20 liquid water content fields with similar characteristics. Incrementally greater x-z plane 
tilt and horizontal stretching were imposed on each of the 20 fields, so that an ensemble of scenes 
was produced for each level of distortion. Broadband solar radiative tranfer was then computed for 
each member of these ensembles using a Monte Carlo transfer model. The domain-averaged solar 
fluxes, diffuse flux ratios, and atmospheric heating rate profiles for the scenes were calculated for 
several combinations of solar zenith and azimuth angle. 

Both tilt and horizontal anisotropy (stretching) were found to significantly affect instantaneous 
domain- averaged solar fluxes, but tilt had more radiative impact than horizontal anisotropy. For 
example, the maximum change in the downward solar flux at the surface caused by tilt was - 
39.9 W m~ 2 ; for stretch it was -12.7 W m~ 2 . Changes in flux caused by cloud anisotropy strongly 
depended on both the solar zenith and azimuth angles in addition to the degree of anisotropy. 
For the stretched scenes, solar transmission increased when the direction of solar incidence was 
parallel to the stretching direction and decreased for broadside incidence. Transmission decreased 
for all solar azimuths except 0° for the tilted scenes. At 0° , transmission increased when the 
clouds tilted towards the sun close to the solar zenith angle. Changes in albedo were the opposite 
of the transmission changes. Although the magnitudes of the changes in the domain- averaged 
absorbed solar fluxes were small (maximum +5.0 W m~ 2 for tilt and +1.3 W m~ 2 for stretching), 
the spatial distortions of the cloud fields caused significant alterations (up to 0.65 K day -1 for tilt 
and 0.25 K day -1 for stretching) in the atmospheric heating rate profiles in many cases. 

Allowing the solar azimuth angle to vary randomly did not eliminate the effects of cumu- 
lus anisotropy, particularly vertical tilt. For random solar azimuths, the maximum change in the 
downward solar flux at the surface caused by tilt was -27.5 W m~ 2 while for stretch it was - 
3.5 W m 2 . For scenes with little tilt or horizontal anisotropy, these values were much smaller 
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(-1 to -5 W m” 2 ). It should be noted, though, that applying completely random azimuths is an 
extreme situation, since solar azimuth varies over a limited range even over extended periods at 
most locations. 

Independent pixel and tilted independent pixel approximation (IPA and TIPA) radiative trans- 
fer calculations were also performed for all of the stochastic scenes. Comparison of domain- 
averaged solar quantities computed using the independent pixel approximation and the full three- 
dimensional Monte Carlo model showed that the IPA was unable to achieve reasonable accuracy 
for the tilted and stretched cumulus fields at most angles of solar incidence. This was because 
the bulk of the radiative effect of the distorted cloud fields was due to changes in their effective 
cloud fraction. (By definition, the IPA is incapable of modelling the incidence of radiation on 
cloud sides and so is only sensitive to nadir cloud fraction.) Errors in the TIPA flux results were 
more predictable, being always negative and on the order of 10 W m 2 . The uniformity of TIPA 
flux errors implies that photon transport perpendicular to the direction of incidence depends only 
on basic cloud structure and properties. This makes TIPA more useful than IPA for approximate 
calculations of radiative transfer through anisotropic clouds. A possible parameterization approach 
for the inclusion of anisotropy in TIPA computations when the full 3D structure of the cloud field 
is unknown was outlined. 

Changes in solar radiative propagation due to cumulus anisotropy may be important in some 
atmospheric applications. Although anisotropy effects may be reduced by relative motion of the 
sun over time or variability of cloud shapes within a large cloud field, we estimate that cumulus 
anisotropy could produce 5-12 Wm~ 2 variations in the downwelling surface solar flux over periods 
of 2.5 hr under optimal conditions. 
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Figure 1: Example LES fields used as input to stochastic cloud field generator. For each scene, 
integrated liquid water path is shown for views in three directions. 
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Figure 2: Horizontal anisotropy (left) and vertical tilt (right) in stochastic cloud scenes as a function 
of the corresponding control parameter. Error bars indicate the standard error of the mean over the 
20-scene ensembles. 
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Figure 3: Example set of stochastic cloud fields with increasing horizontal anisotropy. For each 
scene, integrated liquid water path is shown for views in three directions. 
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Figure 4: Example set of stochastic cloud fields with increasing levels of imposed vertical tilt. For 
each scene, integrated liquid water path is shown for views in three directions. 
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Figure 5: Cloud fraction and average liquid water content as a function of height for input and 
stochastic fields. Solid lines represent ensemble means while dashed lines indicate the mean plus 
and minus one standard deviation, a) Tilted fields, b) Horizontally stretched fields. Note that all of 
the mean value curves are coincident. 
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Figure 6: Definition of solar azimuthal angles for Monte Carlo calculations. Azimuths of 0° and 
180° are parallel to the axis of distortion in all scenes, while 90° is broadside to this axis. Clouds 
tilt towards 0° and away from 1 80° . 


38 


Transmission 


0.8 


0.7 


0.6 


0.5 



o 

o 

0 


0 

Cl 


0.3 


0.2 


0.1 


0.0 


I t- • < - I f T j 



sza = 0 az=0 

sza = 30 az=90 

sza = 60 az = random 


0.08 0.15 0.25 0.35 0.45 0.55 

AP 


0.250 


c 0.225 
o 

"q_ 

o 

(/) 

_D 

< 0.200 


0.1 75 

0.08 0.15 0.25 0.35 0.45 0.55 


~j 


sza = 0 — 

sza = 30 — 

sza = 60 — 


az=0 
az=90 
az = random 


0.500 


0.375 
o 

D 

® 0.250 




-- i---t 


i — i 


00 


0.1 25 


0.000 


sza = 0 az=0 

sza = 30 az=90 

sza = 60 az = random 


0.08 0.15 0.25 0.35 0.45 0.55 


AP AP 

Figure 7: Domain-averaged transmission, reflection, absorption, and diffuse flux ratio for cloud 
scenes with increasing horizontal anisotropy. Error bars indicate the standard error of the mean 
over the 20- scene ensembles. 
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Figure 8: Domain-averaged transmission, reflection, absorption, and diffuse flux ratio for cloud 
scenes with increasing vertical tilt. Error bars indicate the standard error of the mean over the 
20- scene ensembles. 
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Figure 9: Differences between heating rate profiles for tilted and untilted cloud cases. The black 
curve represents the mean heating rate profile for the untilted cases and corresponds to the values 
listed on the x-axis. The remaining curves indicate 10 x the difference between the mean heating 
rates for the tilted and untilted cases. 
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Figure 10: Differences between heating rate profiles for stretched and isotropic cloud cases. The 
black curve represents the mean heating rate profile for the isotropic cases and corresponds to the 
values listed on the x-axis. The remaining curves indicate 10 x the difference between the mean 
heating rates for the stretched and isotropic cases. 
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Figure 11: Differences between domain- averaged transmitted fluxes computed using the inde- 
pendent pixel approximation (IPA) or tilted independent pixel approximation (TIPA) and full 3D 
computations for cloud scenes with increasing levels of horizontal anisotropy and vertical tilt. Top: 
IPA - 3D. Bottom: TIPA - 3D. Error bars indicate the standard error of the mean over the 20-scene 
ensembles. 
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Figure 12: Domain averaged flux results for tilted cloud fields as a function of effective cloud 
fraction (left) and mean cloud optical depth (right) from full 3D computations. Top: Transmission. 
Bottom: Diffuse flux ratio. 
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Table 1: Monte Carlo noise estimates for individual 3D runs. 


Domain- averaged values 

Quantity 

Error 

Downward surface flux 

0.03% 

Upward flux, TOA 

0.13% 

Total absorption 

0.05% 


Absorption profiles 

Quantity 

Error 

Average value in cloud 

1.0% 

Cloudy region maximum 

1.4% 

Clear sky maximum 

0.6% 
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Table 2: Differences between downward surface fluxes computed for stretched and isotropic cloud 
cases [Wm~ 2 ]. 



Anisotropy 

0.08 

0.15 

0.25 

0.35 

0.45 

0.55 

SZA Az. 

Flux 

Downward flux difference 

0 rand 

955.7 

0.1 

-0.7 

-0.4 

0.2 

-1.2 

30 0 

801.0 

1.1 

0.9 

2.0 

3.0 

2.4 

90 

801.3 

-1.3 

-3.0 

-4.4 

-4.7 

-7.8 

rand 

800.9 

-0.2 

-1.1 

-1.2 

-1.1 

-3.0 

60 0 

400.5 

2.7 

3.9 

5.7 

7.6 

9.0 

90 

402.2 

-2.7 

-5.4 

-7.6 

-9.4 

-12.7 

rand 

401.0 

-0.1 

-1.1 

-1.6 

-2.0 

-3.5 
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Table 3: Differences between domain- averaged absorbed fluxes for stretched cloud and isotropic 
cases [Wm~ 2 ]. 



Anisotropy 

0.08 

0.15 

0.25 

0.35 

0.45 

0.55 

SZA Az. 

Abs. 

Absorbed flux difference 

0 rand 

274.4 

-0.1 

0.1 

0.1 

-0.0 

0.1 

30 0 

248.0 

-0.3 

-0.3 

-0.4 

-0.6 

-0.6 

90 

247.8 

0.2 

0.5 

0.8 

0.8 

1.3 

rand 

247.9 

-0.0 

0.2 

0.2 

0.2 

0.4 

60 0 

165.9 

-0.4 

-0.5 

-0.7 

-0.9 

-1.2 

90 

165.6 

0.3 

0.6 

0.8 

0.9 

1.2 

rand 

165.7 

0.0 

0.1 

0.1 

0.2 

0.2 
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Table 4: Differences between downward surface fluxes computed for tilted and untilted cloud cases 
[Wm- 2 ]. 



Tilt 

0.0 

0.33 

0.66 

1.0 

1.33 

1.66 

SZA Az. 

Flux 

Downward flux difference 

0 rand 

955.8 

-2.5 

-8.7 

-17.6 

-26.8 

-36.0 

30 0 

801.0 

5.1 

5.9 

2.7 

-3.3 

-11.3 

90 

801.4 

-1.8 

-6.9 

-14.0 

-21.7 

-29.8 

180 

801.7 

-6.8 

-14.9 

-23.8 

-32.2 

-39.9 

rand 

800.9 

-1.3 

-5.6 

-12.0 

-19.6 

-27.5 

60 0 

400.5 

6.4 

12.2 

16.9 

20.2 

21.1 

90 

402.1 

-0.8 

-2.9 

-6.2 

-10.0 

-14.1 

180 

401.9 

-6.2 

-12.2 

-17.9 

-23.5 

-28.8 

rand 

401.0 

-0.3 

-1.3 

-3.1 

-5.6 

-8.7 
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Table 5: Differences between domain- averaged absorbed fluxes for tilted and untilted cloud cases 
[Wm~ 2 ]. 



Tilt 

0.0 

0.33 

0.66 

1.0 

1.33 

1.66 

SZA Az. 

Abs. 

Absorbed flux difference 

0 rand 

274.4 

0.3 

1.2 

2.5 

3.7 

4.9 

30 0 

248.0 

-1.0 

-1.2 

-1.0 
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3.4 

180 

247.8 
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4.2 
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247.9 
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165.9 
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